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Abstract. In this paper we show that there is a Lipatov bound for the ra- 
dius of convergence for superficially divergent one-particle irreducible Green 
functions in a renormalizable quantum field theory if there is such a bound 
for the superficially convergent ones. The radius of convergence turns out to 
be min{p, l/i>i}, where p is the bound on the convergent ones, the instanton 
radius, and b\ the first coefficient of the /3-function. 



1. The set-up 

1.1. Introduction. In this paper we explore the recursive structure of the short- 
distance sector of a renormalizable quantum field theory. Such theories have a finite 
number of distinct amplitudes r G TZ C A which need renormalization. We decom- 
pose each Green function in accordance with structure functions which need renor- 
malization, the remaining structure functions and superficially convergent Green 
functions contribute to amplitudes a € A \ TZ in the set of all amplitudes A [M] . 

For such a theory, the Dyson-Schwinger equations give any amplitude ,s G A = 
A s TZ U TZ in terms of all those amplitudes r G 7Z in the theory which need 
renormalization, and in terms of an infinite series of integral kernels, the skeleton 
graphs of the theory. The very fact that we can renormalize by modifying the 
Lagrangian implies that we have a sub Hopf algebra at our disposal which has the 
graded elements c r k as generators which correspond to the fc-th order contribution 
to the amplitude r [T51 H"4] . 

We can relate the integral kernels above to the primitives of the Hopf algebra 
structure underlying renormalization. The growth of these primitives is hence deter- 
mined by the growth of integral kernels provided by overall convergent Green func- 
tions. A bound on such a growth is hence obtainable from a bound of amplitudes 
in A \ 1Z, where results in constructive field theory are principally available. The 
typical example is quantum electrodynamics (QED), where the primitives for the 
vertex function are given by the superficially convergent four fermion e + e _ — > e + e~ 
scattering kernel, two-particle irreducible in a suitable channel. 

To proceed from there to a bound for amplitudes in TZ, we need a handle on the 
behaviour of the singular integrations which encompass the short-distance singular 
sector. Here, we proceed by chosing a suitable set of primitives such that we 
can reduce the Dyson-Schwinger equations to recursive equations acting on one- 
variable Mellin transforms. A simple study of the conformal symmetries in the 
corresponding primitives suffice to determine the form of these Mellin transforms. 



Date: December 17, 2006. BU-CMP/06-05. 

*Research supported by grant nsf-dms/0603781 Supported by CNRS. 

1 



2 



DIRK KREIMER AND KAREN YEATS 



For any r € TZ, we can define a Green function 

(1) G r (a,L)^l±Y / ll(a)L k , 

k 

where L = In— Q 2 /[i 2 measures the scaling behaviour under the renormalization 
group flow with respect to to a single Euclidean kinematical variable Q 2 < 0. 

Hence, if we have a recursive set of equations giving Green fuctions G r (a, L) in 
terms of themselves, inserted into integral kernels, we can apply Mellin transforms 
as defined below upon using 

(2) G r (a,ln-k 2 /^) = G r ( a) d p )(~Y 

to reduce the evaluation of all integral kernels above to a study of corresponding 
Mellin transforms. 

In such a situation, we hence show how the growth of superficially divergent 
amplitudes is related to the growth of the superficially convergent ones. To the 
extent the latter is under control by results in constructive field theory, we get 
results for the former. 

The representation of primitive integral kernels through their Mellin transform 
then allows us to turn the Dyson-Schwinger equations into a recursive system which 
determines the 7^ in terms of the 71 , and determines 71 (a) recursively through the 
Taylor coefficients of the relevant Mellin transforms. Modest knowledge of the 
structure of those transforms: 

(3) M(p){p) 



P(l-P)' 

and a Lipatov bound for them, al p 'res„ ~ |p|! C l p l for some suitable c at large \p\, 
allows us to show that solutions to Dyson-Schwinger equations for Green functions 
G"*(a, L) have a similar Lipatov bound. 

We emphasize that our construction of a basis of primitives with a given Mellin 
transform resolves overlapping divergences, thanks to the Hochschild cohomology 
of the relevant Hopf algebras [ID] . 

We have thus a self-similar recursive system determining the formal sums r r , 
r 6 TZ in terms of themselves and the action of suitable maps B^ r = X)| P |=fc -S+ r 
[H]. For all r £11, 

(4) X r (a) =I±J2a k B+ r (X r Q k ) 1 

k>l 

where aQ is the invariant charge of the theory, and Green functions are obtained 
as 

(5) G(a,L) = 4> R (X r (a))(L), 

for renormalized Feynman rules 4>r- 

The study of the Hochschild one-cocycles B+ l ' r is crucial for a QFT. The Hopf 
algebra elements B+ r (T) provide the very integral kernels above underlying the 
DSEs and are the terms which drive the recursion. Their consistent construction 
automatically takes internal symmetries of the theory into account, [13] . by dividing 
the Hopf algebra by suitable Hopf ideals [231 (HI 03] ■ A further ideal is defined by 
our desire to concentrate on amplitudes in TZ. 
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1.2. Ideals. Green functions in field theory decompose into structure functions 
which have logarithmic short-distance singularities. Locality of field theory guar- 
antees that these short distance singularities are invariant under changes of dimen- 
sionfull parameters in the theory, allowing for local counterterms. This is not the 
only invariance which we can observe in counterterms. A chosen subgraph needs 
the same counterterm wherever it appear in a larger graph, and in whatever ori- 
entation it is inserted, as long as we work in a symmetric renormalization scheme. 
Furthermore, its counterterm is invariant under modification of its external mo- 
menta. We can hence isolate short-distance singularities in subgraphs in a manner 
such that all subdivergences are functions of a chosen internal momentum of the 
cograph under consideration. 

In so doing we enlarge the Hopf algebra of graphs into a Hopf algebra of colored 
graphs or decorated rooted trees, which makes the underlying Hochschild cohomol- 
ogy obvious and resolves overlapping divergences [TU1 H2] ■ 

We emphasize that the choices above depend on the theory, renormalization 
scheme and physical problem one wants to study. Here, we are only concerned with 
the general set-up, and will only sporadically restrict to a specific theory. Our final 
result applies to any renormalizable field theory. Specifics of any particular such 
theory are to be discussed elsewhere. 

The freedom we have in the above choices can be summarized by saying we work 
in an ideal given by 

(6) 5^r // o i End _1 r / = 0, 

i 

with 

(7) EndO^) = End(r!)End(r 2 ). 

Here, we use Sweedler notation A(r) = Yli F' <8> T" and the fact that Feynman 
graphs have an operadic structure, so that for each term in the sum we have an 
operadic composition [12) of graphs such that 

(8) r = r" o, r'. 

Finally, End is a map which implements such a choice: it chooses a configuration 
of external momenta for the subgraph T', permutes the orientation of external legs, 
or changes the insertion place: any modification which will not modify the required 
counterterm is allowed. By construction, elements in our ideal are superficially 
convergent, and hence contribute nothing to the short distance singular sector. The 
above ideal is determined by the choice of the map End, which is then determined 
by the physics one wants to study. 
We then decompose graphs as 

(9) r = J2 r " °i End- 1 (r')P(r") + J2 r " °* End- 1 r'F(r") J p(r'). 

i i 



By construction, the first term in the rhs is a primitive element p = p(T) in the 
Hopf algebra of decorated rooted trees, and End -1 = Endo S. 

Note that for a primitive graph 7 and an endomorphism End which changes its 
external leg structure to say zero momentum transfer, the above ideal simply says 
that the short distance singular sector of 7 and End(7) are the same. 
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More instructive is the example of the case of say a graph with a single divergent 
subgraph r = T" o V . We then have 

(io) P (T) = r" o [r' - End(r')] . 

This is clearly primitive as [r' — End(r')] is in the chosen ideal, and p(T) in ([9]) is 
primitive in general because all its subgraph are in the ideal, by construction. We 
emphasize that in all our applications we will never discard any terms on the rhs 
of Q, but will always calculate the full graph T, but use jl| as a convenient tool 
to come to a manageable form for the primitives which drive the recursion, on the 
expense to have an enlarged set of such integral kernels. 

For our purposes, we choose an ideal which sets to = in structure functions and 
evaluates external momenta at a symmetric point — Q 2 with regard to their external 
momenta, and iterates graphs into each other in accordance with the definition 
of the Mellin transform defined below. As a renormalization scheme we choose 
subtraction at a fixed — Q 2 = /i 2 . 



1.3. Mellin transform. Any so constructed primitive p = p(T) is a degree ho- 
mogenous combination of Hopf algebra elements of degree \p\. We let — Q 2 be the 
above kinematical variable which we keep. We define the Mellin transform M (p) 
of p as 



(11) M(p)(p) 



[-Q 2 Y J Int p (-Q 2 ) | -1 f>?r I fl d 4 k t , 



where Int p (— Q 2 ) is the integrand determined by p. We let 

(12) Int;(-Q 2 ) := Int^-Q 2 ) - Intp^ 2 ). 

We define renormalized Feynman rules for a symmetric momentum scheme with 
subtractions at — Q 2 = /j, 2 by 




(13) MB p + (X))(-Q'/^) = Int-(-Q^) I > 6 R {X){-ki/^ 



We have 

(14) MBl(X)){-Q 2 1 > 2 ) = lim MX){d- p )M(j>){p) [[-Q 2 /^T P l] , 
where <9_ p = 

A Green function is then defined as the image under such Feynman rules applied 
to a fixpoint of a combinatorial DSE. In the coming sections, we first study the 
asymptotics of a combinatorial DSE, then the growth after applying <f>R. 

Let us now look at an example in (f)\ theory. We define the vertex Green function 
by setting —Q 2 = s = t = u, and set m = in all propagators in accordance with 
an ideal which isolates the short distance singularities in massless Green functions. 

Taking the symmetry in external legs into account, the Hopf algebra series reads 
to order g 2 

as) 1+4 o +*9 2 (\ oo +\ 
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We define a primitive 

P2 = \ OO + 

(is) = \ OO - 



i r 

2 
1 

2 



O °End( O 



We define a combinatorial DSE 
(17) X(g) = l + g^(X(g) 2 ) 



where p\ = (^} ■ 

It has a fixpoint which reproduces to order g 2 the above expansion. 
Let us turn this into an integral equation: 



G{gM-Q 2 /p 2 ) 



G(g,\n-k 2 /p 2 ) 2 



d 4 k 



(18) 



3 5 2 



1 



k 2 (k + Q) 2 
[G(g,]n-k 2 /» 2 )+G(g, In 
2k ■ I - 2(k + I) ■ s - 2s 2 



W)] 3 

d 4 kd 4 l, 



k 2 l 2 (k + Q) 2 (l + Q) 2 (l-k) 2 

where }_ indicates subtraction at — Q 2 = p, 2 . One verifies that the solution to 
this integral equation agrees to order g 2 with the perturbative renormalization in 
a symmetric momentum scheme for the vertex function in <j>\ theory, as it should. 
Its solution also agrees to any order in g 2 with the leading order in L = In — Q 2 /p 2 
in that order of g 2 , and also with the next to leading order, as we have taken the 
two loop primitive into account. 

Addition of further primitives delivers the lower powers in L. The Mellin trans- 
form Mr,, = —rr — 7 
Pi p(i-p) 



is trivial, the Mellin transform M P2 much less so and can be 



found in [2]. 

1.4. Mellin transforms as a geometric series. The above Mellin transforms 
M p (p) have the form 

Up) 



(19) 



M p (p) 



P(l-P)' 

in a natural manner, where the pole at p = reflects the short-distance singularity, 
and the pole at p = 1 reflects the fact that in our ideal we have massless internal 
propagation. 

The series f p (p) — res p + O(p) determines the residue res p of the transform. 
We keep the factor 1/(1 — p) explicit though, to maintain the form of the Mellin 
transform similar, in accordance with its conformal symmetries at —Q 2 = 1. 

By the above definition of the Feynman rules it is a purely algebraic exercise 
[21 [TT] to define a new series of primitives pi — p, p2 — B+(B+(T)) 
and so on, such that 



M p(p) =^2<ln 



rcs„ 



(20) ^ pvr , 

We henceforth assume a basis of primitives {p} such that 

(21) M p (p) 



res,. 



P(l-P)' 
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which is convenient in the following. This finishes the discussion of our analytic 
set-up and we next remind ourselves of some basic properties of recursive systems 
relevant to Dyson-Schwinger equations. 

2. The Universal law 

2.1. The classical case. A classical result of combinatorics is Polya's asymptotic 
formula for the number, t(n), of rooted, unlabelled trees with n vertices |20j (trans- 
lated in [21]): 

(22) t(n) ~ Cp- n n- 3 / 2 

where C is an explicit constant and p is the radius of convergence of the generating 
function T(x) = J2 n >i t(n)x n 01 the class of all rooted, unlabelled trees. Thus p is 
the reciprocal of Otter's tree constant [19] (sequence A051491 in [22]). The key to 
Polya's proof is to convert the recursive equation 

T(x) = xexp ( T(x m )/m j 

^m>l ' 

for rooted, unlabelled trees to a bivariate function 

E(x, y) = xe v exp ( T(x m )/mj 

^ m>2 ' 

at which point the recursive equation becomes T(x) — E(x, T(x)). Since p < 1 this 
rewriting uses the new variable y to isolate the portion of the recursive equation 
which controls the radius of convergence. 

Consequently, y — E(x, y) is amenable to Weierstrass preparation, which tells us 
that around (p, T(p)), y — E(x,y) is a product of a holomorphic function which 
is nonzero around (p,T(p)) and a monic polynomial in y with coefficients which 
are analytic in x. The degree of the polynomial is the order of the first nonzero 
derivative. The failure of the implicit function theorem at p gives that 

T(p)-E(p,T(p))=0 

l-(d y E)(p,T(p))=Q 

-{d yy E){p,T{p))£Q. 

So the Weierstrass polynomial is quadratic which gives that T has a square root 
singularity at p: 

T{x) = f(x)+g(x)yfp=Z 
with / and g analytic at p. Then the Cauchy integral theorem gives the desired 
asymptotics along with a formula for C . 

Asymptotics of the form C p~ n n~ z / 2 have since been widely found for classes of 
rooted trees with recursive definitions and hence (|22p is sometimes known as the 
universal law. Some notable examples generalizing the reach of Polya's analysis in- 
clude [3], 0) and [T7]- Surveys include [5] and [TB]- See also chapter VII of Flajolet 
and Sedgewick's forthcoming book Analytic Combinatorics currently available as 
an online draft [70 

For all but the last of the above the focus is on improved analytic conditions 
which still must be checked for any particular case of interest. One approach for 
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further generalization is to specify certain combinatorial constructions and then use 
these constructions to build recursive equations with solutions which automatically 
satisfy the universal law. This is the approach of pQ. Flajolet and Sedgewick also 
use a framework of combinatorial constructions to build their recursive equations, 
getting the universal law for various schemas, such as the exp-log schema [7J VII. 1]. 

2.2. Recursive systems. Another approach to generalizing Polya's analysis is to 
move to recursive systems of equations while restricting the complexity of each equa- 
tion. This is of particular relevance to applications to counting Feynman diagrams 
since generally only polynomials and geometric series are needed, but systems are 
difficult to avoid. 

First note that any geometric series can be converted to a polynomial at the 
expense of a new variable. Namely replace 1/(1 — X) with a new variable F and 
add the equation 

F = 1 + F ■ X. 

In view of this reduction we can focus on polynomial systems. 

Nonnegative polynomial recursive systems give the universal law under some 
reasonable conditions, as was shown independently by Drmota 6J, Lallev[16]. and 
Woods [25 . For our purposes the full generality of the above are not necessary and 
we'll follow the presentation of Flajolet and Sedgewick [7, VII. 6. 3]. 

Suppose 

yi = <^ 1 (x,y 1 , . . .,y m ) 



with the $i polynomials with real coefficients. 

There are five conditions which together guarantee that each component solution 
to this system satisfies the universal law. First we say the system is nonlinear if at 
least one of the <frj is nonlinear in y\, . . . , y m . The system is nonnegative if each $^ 
has nonnegative coefficients. 

The next condition guarantees that the system does in fact behave recursively. 
For y = (yi, ■ ■ ■ ,y m ) € R^a;]]" 1 define the x- valuation by val(y) = min^val^)) 
where the valuation of a series picks out the index of the first nonzero coefficient, 
i.e. yal(Y] n _ k a n x n ) = k when ^ 0, and with the convention that val(0) = oo. 
Also define d(y,y') — 2~ val ( 2/ ~^ ). Then the system is proper if 

d($(y), $(y')) < Kd(y,y') for some K < 1. 

To guarantee that the system behaves as one system rather than many we need 
to guarantee that all variables play a role in each equation. Specifically, define the 
dependency graph of the system to be the directed graph on {1, . . . , to} with an 
edge from k to j if yj figures in a monomial of y). The system is irreducible 

if the dependency graph is strongly connected, that is there is a path between any 
two vertices following edges only in their forward direction. 

Finally, to remove spurious zero coefficients in the solution series, and to avoid 
extra singularities on the circle of convergence, define a power series T(x) to be 
aperiodic if it is not the case that there is a power series U(cc) and integers a > 
and d > 2 such that T(x) = x a XJ(x d ). If the solutions for each yi of the system are 
aperiodic then the system itself is said to be aperiodic. 
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Then 



Theorem 1 ([JJ Theorem VII. 6). Suppose y = <fr(Tj) is a polynomial system that is 
nonlinear, proper, nonnegative, and irreducible. Then all component solutions yj 
have the same radius of convergence p < oo and have a square root singularity at 
p. If furthermore the system is aperiodic then all yj satisfy the universal law. 

An outline of the proof is as follows. For more details see [7J Theorem VII. 6]. 
Properness gives that the component solutions are unique, because the solutions can 
be generated recursively from the zero vector, and along with nonnegativity we get 
further that the component solutions have nonnegative coefficients. Irreducibility, 
along with nonnegativity to prevent cancellations, gives that the radius of each yi is 
the same value p. Aperiodicity forces p to be the only singularity of each component 
solution on the circle of convergence. 

The key to the proof is the square root singularity at p. Each of the three in- 
dependent proofs takes a different approach, though in all cases, as in the single 
equation situation, it comes down to using the failure of the implicit function the- 
orem at p. Nonlinearity is necessary to avoid singularities which themselves are 
linear rather than of square root type. 

Drmota [6] proceeds by solving a subset of the equations for the remaining vari- 
ables and then substituting back into the remaining equations to reduce the number 
of equations in the system. Iteratively he is able to reduce to one equation at which 
point the system can be treated classically. 

Lalley [16] , summarized in [JJ Theorem VII. 6], considers the linearized Jacobian 



and uses Perron-Frobenius theory to show that the largest eigenvalue in absolute 
value at x = p of the Jacobian is precisely 1 and that there is an eigenvector v with 
positive coefficients. Then multiplying the system by v and expanding around p 
gives the desired asymptotics. 

Woods [25] also uses the Jacobian and the largest eigenvalue 1. He continues 
the analysis on block upper triangular matrices in order to deal with certain non- 
irreducible cases. 



3.1. QED with one primitive per loop order. Massless quantum electrody- 
namics provides three monomials in the Lagrangian which need renormalization, 
corresponding as one-particle irreducible Green functions to the inverse fermion and 
photon propagators, and the vertex. We do not yet impose the Ward identity and 
proceed by setting up the combinatorial structure as it rather typically holds for a 
theory with three-valent vertices and two types of edges (the degenerate case of a 
single type of edge is easily interfered). 

The smallest Hopf algebra which allows for a correct renormalization of the full 
theory treats the sum of all primitive graphs at a given loop order as a single 
Hochschild cocycle, and we start from there. This situation is described by the 




3. The universal law for Feynman diagrams 
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following system. 



X 2 = x 
X 3 = x 



2k+l 



^~ (i-x 2 )^(i-x 3 y 



(1-X 2 )(1-X 3 ) 



[l-x 2 f 

After converting the geometric series we get 



'x, 


= 1 + X1F1 


Fx 


= xX\FlF 3 + xX\FlF 3 F 


x 2 


= xX 1 F 2 F 3 


F 2 


= 1 + F 2 X 2 


X 3 


= xX x F$ 




= 1 + ^3*3 



<£> is nonnegative, polynomial, nonlinear, and irreducible. $ can be seen to be 
aperiodic simply by calculating the first few terms. $ itself is not proper. However 
<1> 2 describes the same solutions and the other properties, nonnegative, polynomial, 
nonlinear, irreducible, and aperiodic, remain true for powers. 
To see that <£> 2 is proper suppose we have 2 vectors 



v = (xi, fl,X 2 , f 2l x 3 , f 3 ) 



and 



v — ( x it fii x 2i /21 x 3> /a) 
at distance 2~"; so write x\ = x\ +x'[ with x'[ having no term of degree less than n 
and similarly for the other coordinates. Consider the difference in X\ coordinates 
after applying $: 

r I pi II nil II p pll 

xifi - x 1 f 1 = -x 1 f 1 -X 1 fi- XlJx 

which has no terms of degree less than n. Further if /1 has no constant term and 
/{' has no term of degree n then the difference in the X± coordinate has no terms 
of degree less than n + 1. Argue similarly for F 2 and F 3 . For the X 2 coordinate 
after applying $ we get a difference of: 

W2/3-ZZ1/2/3 = -^1/2/3 -xxZfah-xxififs-zxifSfa '-^1/2/3-^1/2/3' 

which has no terms of degree less than n+1. Notice also that the new X 2 coordinate 
has no constant term for both initial vectors. Argue similarly for A3 and F\. 

Now consider applying $ a second time. Apply the above arguments again but 
notice that we are in the "further" case for X\, F 2 and F 3 , so all coordinates now 
have a difference with no terms of degree less than n+1, that is the distance has 
decreased by at least 1/2. Thus any 1 > K > 1/2 will give that $ 2 is proper. 

Consequently by Theorem Q] we know that all 6 solution series to <I> 2 and hence 
all 3 solution series to the original system have the same radius of convergence p 
and have coefficients with the asymptotic form 



Cp 



-n n -3/2 
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In order to understand the asymptotic growth rate for the QED system with 
primitives summed at each loop order it remains to understand the radius p. Sub- 
stituting X2 and X3 into F2 and F3 respectively we get F2 = 1 + xXiF^F^ = F3. 
Thus the original system can be rewritten 

xXfFi 



Xl ~ 1 ' 1 - xXlFl 
F 2 = l+xX 1 Fl 



Rearrange 
then solve to get 



A 1 = l 



F 2 = 



l-Xi(F 2 -l) 



1 - 2Xl 



X 1 (1-2X 1 ) 

Substitute back into the equation for F2 and expand to get 

-x + Xi + (6x - 5)X% + 8X'f + (-12a; - i)Xf + SxXf = 
As a polynomial in X\ this has discriminant 

4096a; 2 (32a; 2 - 8a; + l)(-2 + 27a;) 2 
So the radius of the system is 

2 

27 

In view of the fact that the radius is an important value associated to the system 
and that this system is canonically associated to QED we're led to the following 
question as to what is the physical meaning of 2/27 in QED? The relevance of 
this number extends to any system which provides a recursive system similar to 
QED. We can proceed similarly for other theories; the results of some examples are 
summarized in Appendix [A] 

3.2. Polynomially many primitives per loop order. Combining all primitives 
at a given loop order into one Hochschild cocycle which drives the Dyson-Schwinger 
equation defines the smallest sub Hopf algebra which still renormalizes the full 
theory correctly. It is often instructive to disentangle the primitives in different 
ways, for example in accordance with the transcendental nature of res p . This 
motivates to consider a slightly more general condition on the number of primitives, 
enough to apply the polynomial systems result, we next assume there to be p(k) 
primitives at k loops where p is a polynomial. 

To see that this circumstance reduces to a nonnegative polynomial system it 
suffices to show that 

k>e 

can be written as a sum of powers of geometric series. This follows from two facts. 
First, the falling factorials 

{k{k- l)---(k- n+ 1) :n>0} 

form a basis for polynomials in fc, and 

n 

fc" = £s 2 (n,j)fc(k-l)---(fc-i + l) 
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where S^tn, j) are the Stirling numbers of the second kind, A008277 in [22 , so in 
particular are nonnegative, and hence for p{k) with nonnegative coefficients we only 
need nonnegative coefficients of the falling factorials. 
Second, notice that for £ > n 

^2k(k-l)---(k-n+l)B k 

k>t 

Jn 

= B nJ_y B k 

dB n ^ 



n 

3=0 



= E^-D-(^-"+i + l)Q)(l^Bj 

where all the coefficients are nonnegative. 

One case where there is a natural interpretation is QED with a linear number of 
generators, namely 

^1 = 1 + V#)xS rkn rr 

with Xi and A3 as before and with p(k) linear, which corresponds to counting with 
Cvitanovic's gauge invariant sectors 0]. 

3.3. Other systems. Johnson, Baker, and Willey [9] use gauge invariance to re- 
duce the QED system to 

/ _ \ k 



fc>0 



x(l-X) 
-XI l-X-x 



While amenable to the universal law analysis, this recursive equation can be solved 
exactly by the quadratic formula. We get 



x _ 1 + y 7 ! - ^ 
2 

giving the Catalan numbers, A000108 in [32] , as coefficients. The radius is 1/4 which 
is considerably larger than 2/27, showing how powerful gauge invariance is. Note 
that it is only the inverse photon propagator 1 — X which needs renormalization, 
and that it appears in the denominator. 

4. The growth of 71 
After these considerations of the combinatorial side, we discuss analytic aspects. 
4.1. The recursions. Consider the Dyson-Schwinger equation 



x(x) = 1 - E E * k pmB k +{xQ k ) 



k>l »=0 



where Q — X r with r < an integer, and Pi(k) coefficients, not necessarily polyno- 
mial. We'll use the notation Fk : i{p) for the Mellin transform of the integral kernel. 
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Starting with r < is justified in light of the last example of the previous section, 
we will generalize this soon enough. 

Specializing [T51 (26)] to this case we get the recursion 

(23) 7 fc (x) = -7i(x)(l + rxd x )~f k -i(x) 

k 

independently of the Pi(k). Using the dot notation, j-U — lkU k , of [15] we have 

(24) 7 -i = EE x k Pi(k)(l + 7 • d- p )- rk+1 (l - e- L »)F kii (p)\ p=0 

k i 

Taking one L derivative and setting L to we get 

(25) 7l =^^ 2 ;V i (fc)(l + 7-5- P )- rfe+ V^M(p)lp=o 

k i 

Restricting to pFk t i{p) — — p) allows us to write a much tidier recursion 

for 71. Taking two L derivatives of (f24|) and setting L — we get 

2 72 = - E E ^(fc)(l + 7 • 9-,)- rfc+ yi^(p)U=o 

k i 

= + 7 • d- P r rfe+ v^(p)ip=o + 5>v,*K(fc) 

k i k 

= -71 + E E r fe,^( fc ) a;A; from (ESD 

fc>l i 

Thus we can from now on ignore the sum over i and let p(k) — rk y iPi{k). Then 
from (USD 

7i ^X!^ -2 ^ = Ep( fc ) xfc - 71(1+^)71 

/c>i /c>i 

giving 

Proposition 2. 

n-l 

7i,n =K n ) + E(~ r -? _ 1 )7lj7l,«-J 
i=i 

4.2. Finding the radius. We see from the proposition that if J2p(k)x k is Gevrey- 
n, that is x k p(k)/(k\) n converges, but not Gevrey-(n — 1), then 71 is at best 
Gevrey-n. 

Of most interest for our applications are the cases where only finitely many p{k) 
are nonzero but all are nonnegative and where p(k) — c k kl giving the Lipatov 
bound; so ^p(k)x k is Gevrey-1. Assume that p(k) > and 

k>l 

has radius < p < 00 and f(x) > for \x\ < p. The above two cases are included 
as f{x) a polynomial and f(x) = cx/(l — cx) respectively. In such circumstances 
7 i is also Gevrey-1 and the radius is the minimum of p and — l/(rai) (where we 
view —l/{ra\) as +00 in the case a\ = 0) which we can see as follow^ 



2 Similarly ^2 p(k)x k Gevrey-n for n > 1 leads to 7„ Gevrey-n 
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Let a n = Ji, n /n\. Then ai = 71,1 =p(l) and 

n-l / \ -1 



— , 



/ \ _ 

Pi " ' , 1 ■ , i ■ \ Y\f n \ 

2 Z^(^ r -? - 1 - r(n - .7) - 1) J ajan-j 



= T + h- 1 )E i ^ 

i=i VJ/ 

To achieve an upper bound on the radius of convergence of a nX n take the first 
and last terms of the sum to get 

p(ri) n-2 

a n > — : r oio„_i 

n\ n 

for n > 2. So the radius of ^ a„x" is no more than the radius of the recursively 
defined series with equality above, say 

p(n) n-2 
o n = — ; r Oib„_i 

for n > 2 with &i = a\. Immediately we see that if a\ = the radius of ^ b(n)x n 
is p. Otherwise consider 

„(„ _ i) bn = p(n)n( " - 1} - r (n - l)(n - 2)ft 1 6„_ 1 
n! 

Equivalently with B(x) = b(n)x n we get 

B"(ar) = /"(a) - r&ixB"(a;) 

Solving for B"(x) 

1 + raia; 

which, since differentiation does not change the radius of a series, has radius 
min{p, — l/(ra\)}, and thus so does B(x). 

For the lower bound on the radius we need a few preliminary results. First a 
simple combinatorial fact. 



1 fn\ 6-3+ 1 
> 



Lemma 3. Given < 6 < 1 

n \j 

for 1 < j < On and n>2. 

Proof. Fix n. Write j = An, < A < 6. Then 



1 /n\ _ 1 / n \ n Xn - 1 _ X~ Xn+1 

n \j J n\XnJ ~ (Xn) Xn Xn ~ j 



a 



Second we need to understand the behaviour of a nX n at the radius of conver- 
gence. 

Lemma 4. Using notation as above and with A(x) = a n,x n with radius of con- 
vergence p a we have that limsup a ._ >po A(x)(l + xra\)/ f{x) < 1 
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Proof. Take any < 9 < 1/2. Using the previous lemma 

n— 1 / N -1 



_p(n) 

tin - 



/ \ 



l<7<0n VJ/ $n<j<n-6n yJ 

p(n) ^ ._ 1 /"V 1 " 

^— ~ r i° ajan - j ~ r {eJ 2 ^ aj ' a "^ 

l<j<6n V ' en<j<n-6n 

l<j<8n 8n<j<n-8n 

so the coefficients of A(x) are bounded above by the coefficients of 

f(x) - xrA'(6x)A(x) - £ ( x^A 2 ) (0 e x) 
2 \ ax ) 

Since all coefficients are nonnegative, for any < x < p a we have 

A(x) < f(x) - xrA'(6x)A(x) - T - (x^-A 2 ) (6 e x) 

which is continuous in 8, so for fixed < x < p a we can let 9 — > giving 

A(x) < f(x) — xraiA(x) 

so 

A(x){\+xra 1 ) 
fix) 

for < x < p a - The result follows. □ 

Let p a be the radius of a n x n . If p a = p then we're done, so suppose p a < p- 
To get the lower bound on p a it remains then to prove that p a > — l/(ra\) when 
ai 7^ and to prove a contradiction in the case that a\ = 0. Take any e > 0. Then 
there exists aniV>0 such that for n > N 

n-2 



p(n) 1 X- 

a n < — : raia„_i - r > a^a n - 

n\ n - 1 *r~L 



J=2 

n-2 



n! ^ 



Define 
a 



j=2 

< — ; ra\a n -\ + e } aja n - 

p(") 



■r ^ P\ n ) i v-^ n ~ 1 

ii a n > - rc!C„_! + c^c^ 
l^nT — rc i c n-i + e Sj=i c j c n-j otherwise (in particular when n > N) 

In particular ci = oi. The radius of ^ a„x™ is at least as large as the radius of 
C(x) = °nX n . Rewriting with generating series 

C(x) = f(x) - r ai xC(x) + eC 2 (x) + P e (x) 
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where Pe{x) is some polynomial. This equation can be solved by the quadratic 
formula. The discriminant is 

(26) (ra lX + l) 2 - 4e(f (x) + P e (x)) 

By assumption we are interested in |x| < p where (|26l) has no singularities, so the 
radius of C(ir) is the closest root to of (|26| ; call it p e . Consider e = giving 

(raix + if 

which has — X/(rai) as its closest root to when a% ^ 0; call this value po- 

Then to get the lower bound on the radius of ^2 a n x n it remains only to prove 
the following lemma. 

Lemma 5. With notation and assumptions as above if ai ^ 0, lim e ^o \Pe\ > Po, 
while if ai = we have a contradiction. 

Proof. In view of Lemma [4] this is a short exercise in analysis. 

By construction the coefficient of x n in P e (x) is bounded by a n + rc\c n ^\ < 
a n + ra\a n -i since c n > a n for all n > 1. So P t (x) has coefficients which are 
nonnegative and bounded by those of A(x)(l + rxai). Thus by Lemma |U the 
continuity of P e (x) at p a , and the assumption that p a < p, we see that f{p a ) + 
Pe{Pa) < f(Pa) + lim inf :c ^p a (A(x) (1 + rxai)) < oo. By the nonnegativity of the 
coefficients of / and P t we can choose M > such that \f(x) + P e (x)\ < M 
independently of e for \x\ < p a - 

Suppose a% ^ 0. Take any r\ > 0. Consider \x\ < p a . Choose 5 > such that 
(raix + l) 2 < S implies \x - p \ < r\. Pick e < S/(4M). Then {ra\p f + l) 2 = 
4e(/(p e ) + Pe(pe)) <Sso \p t ~p a \<v 

Suppose on the other hand that a\ = 0. Take < 5 < I . Then, since \p e \ < p ai 
we get that for e < 5 /{AM), 1 = ie(f(p e ) + P e (p e )) < S which is a contradiction. □ 

Taking the two bounds together we get the final result 

Theorem 6. Assume ^2)~>iX k p(k)/k\ has radius p. Then Y]x n ^i. n /nl converges 
with radius of convergence min{,o, — I/(r7i.i)}, where — l/(? , 7i i i) is interpreted to 
mean +oo in the case 71,1 = 0. 

4.3. Nonnegative systems. Now suppose we have a system of Dyson- Schwinger 
equations 

k>l i=0 

for r 6 1Z with TZ a finite set, pl(k) > 0, and where 

q = J] x r ( x yr 

with integers s r < for all r € TZ. 

Then as before from p75|. (26)] we have 

(27) rt(x) = i ^(^) + E s n i(x)xd}\ f k _ x {x) 

again independent of the p r (k). 

Assume that there is one insertion place, and so one variable p, and that the 
Mcllin transform of the integral kernel is a geometric series pFJI i (p) = r^ i:r / (I — p). 
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Rewriting the system of Dyson-Schwinger equations with dot notation we have 
(28) 7 r • i = E E (*0 II C 1 + ^ ' d- P r s > k+1 (l e-^FZMl^ 

k i jell 

As before we can find tidier recursions for the "/{ by comparing the first and second 
L derivatives of ([28]), We get 



7i 

and 



= E E xk pi( k ) II c 1 + ^' • 3-pr 3fe+ VL(p)ip=o 



2 ^ = - E E ^tf (*) n c 1 + ^ ■ d - P )' s ' k+1 p 2F Up)\p-o 

k i jGTZ 



= ~7i + E E r k^rP r l {k)x k since pF kii (p) = — 

k>l i P 

Thus letting p(k) = J2i r k,i-rPi{k) and using (J2TJ) 

k>l k>l jeK 

giving 

n — 1 n— 1 

7i,„ = P r («) + E( -Sri ~ 1 )Ti,iTi,n-i + E E(~ S ^)^1»-^M 

t=l jeK i=l 

To attack the growth of the j{ we will again assume that 

fe>i 

has radius < p r < oo and / r (a:) > for |x| < p r . We'll proceed by similar bounds 
to before. 

Let a r n — r )\ n jn\. Then 



r , v n— 1 / v -1 n— 1 / v -1 

' ,; ' ! ■ E(-^-D«-4 •) +EE(^K-xr) 



rt! 

i=i x 7 jeR «=i 

Taking the last term in each sum we have 

r ^ p r ( n ) ir j 1 71 ~ 2 



r ^ P V' L \ " j 



-a., 



»! i — ' " * ; n 

Let 6^ be the series defined by b\ — a\ and equality in the above recursion. Let 
^ r ( x ) = J2^n xn - Then as before if X^eTC s j a i = tne racu us of B(a;) is p r and 
otherwise consider 

W{x)" = fix)" - (E J ^(x)" 
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Solving for W(x)" we get that the radius of a n xn 1S a t most 

f -1 1 



mm < p r , 



E 



jev 



again interpreting the second possibility to be oo when EjeK s i°i = ^- 

In the other direction take any e > then there exists an N > such that for 
n > N we get 



<<^-[£*xKx+ e E£ 



r 7 

°X-i 



vierc / i=i jerc 



Taking C r (x) to be the series whose coefficients satisfy the above recursion with 
equality for when this gives a result > a r n and equal to a r n otherwise we get 

c r \x) =r(x)- s J a i ] xcr ( x ) + e E cr ( x ) cj (x) + pi{x) 

\i£K / jeV 

where PJ is a polynomial. 

Summing over r we get a recursive equation for J2ren C r (x) of the same form 
as in the single equation case. Note that since each C r is a series with nonnegative 
coefficients there can be no cancellation of singularities and hence the radius of 
convergence of each C r is at least that of the sum. The equivalent of Lemma [4] for 
this case follows from 

E < ± E - E »A E + E E(-^K-*< (" 



^E^-E^iE^ 

rev jev rev 

<-^E<")"YE<-il (e< 



max( 
j - - z — / 10 

i=2 x ' \reV J \reV 



E^r-E s ^iE a - i 

r£R ' 3"£R r€R 



1 

71 1 



+ max(-s> E ( J E °«-< E °i J 

2<i<9n ^ ' \r£R / \reV / 

+ max(- s ,)^ E (") (E a n-)(E<) 

9n<i<n-6n V 7 VrGTC / VrGK / 

for as in Lemma|3]with EreK A r (x) in place of A(x), where A r (x) = ^ a r {n)x n . 
Then continue the argument as in Lemma|3]with EreTC f r ( x ) m pl ace of f{x) and 
maxj(— Sj) in place of — r, and using the second term to get the correct linear part 
as 9 -> 0. 

Thus by the analysis of the single equation case we get a lower bound on the 
radius of a n x ™ °f mm se7j{Ps) — 1/ J2jev s j a {}- I n particular if r € 1Z is such that 
p r is minimal we see that the radius of a r n x n is exactly m.va{p r , — 1/ Ejgtc s i a i}- 
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Suppose the radius of ^ a s n x n was strictly greater than that of ^a r n x n . Then 
we can find [5 > S > such that 

< > p n > 6 n > < 

for n sufficiently large. Pick a k > 1 such that a s k > 0. Then 

5 n > > > „ r ...^ 



s r k\a s k s r k\a s k , 

n---(n-fc + l) n---(n- fc + 1)' 



so 



5\"- fe fc! 
> 



— s r a^ \/3/ n---(n — fc + 1) 

which is false for n sufficiently large, giving a contradiction. Thus all the ^ a s n x n 
have the same radius min re 7j{p r , — 1/ X^gtc s i a i} 

From this we can conclude that each ^2x n j{ n /nl also converges with radius 
mm re fc{p r , — 1/ X^jeTC s j7i ,i}i where the second possibility is interpreted as oo 
when £ jeK s 3 7^ 4 = 0. 

4.4. Systems with some s r > 0. Let us relax the restriction that s r < and that 
p r (n) > 0. It is now difficult to make general statements concerning the radius of 
convergence of the ^ d^x n . For example consider the system 

1/ n n—l / \ — 1 n— 1 / \ 

i_P W , \^, , -n„i„i f^M Y^„-„i„2 f n 



a 



IV J. / \ —l ft -I / \ - 

+ E ( 2 J - ^A-j Q - E j a K -i Q 

n— 1 / \ — 1 n—l / \ - 

X>' + tyX-i Q + E ^X-,- Q 



so si = — 2 and s 2 = 1. Suppose also that 

p 2 (2)-0 
a} = 



p 2 {n) = -2{n-l)\a\al 



-l 



Then a| = and inductively = for n > 2 so the system degenerates to 




a\ if n = 1 
otherwise 



We still have a free choice of and hence control of the radius of the a 1 series. 

On the other hand the a 2 series trivially has infinite radius of convergence. 

Generally, finding a lower bound on the radii of the solution series, remains 
approachable by the preceeding methods while control of the radii from above is 
no longer apparent. 
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Precisely, for any e > 



k: I < 



- E s J a i 

n-2 



n-2 



dl + SK-ari-iJIKIK- 



< 



+ EE 

Ip»I 



E s J' a i 

jell 



n-1 

-il + ^EEKIK-l 



So, for a lower bound on the radius we may proceed as in the nonnegative case using 
the absolute value of the coefficients and achieving that the radius of the in 1 ™! 
and hence that of £ n 7i n / n ^ 1S a ^ least 



mm 

relZ 



Pr, 



s j7i,i 



0. 



where the second possibility is interpreted as oo when J2jen s j'7i 

Note that this gives the lower bound on the radius of convergence as the minimum 
of the first instanton singularity (which one expects to be the radius for p(k)/k\) 
and the inverse of the first term in the j3 function of the theory. Furthermore, we 
emphasize that Ward identities typically allow a restriction to systems where all 
s r < 0. A more detailed discussion will be given in future work where the general 
approach described here will be discussed with regard to the specific details of the 
relevant renormalizable theories of interest. Finally, we note that the apperance of 
the inverse of the first term in the /3-function makes sense: in the conformal case 
of a vanishing /3-function we would not expect a constraint on the minimum of the 
radius to come from perturbation theory. 

5. Applications of the growth of ji 

While expectations for the growth of p{k)/k\ in terms of instanton singularities 
are routine in the context of path integral estimates, the path integral is merely a 
successful heuristic to parameterize our lack of understanding of quantum field the- 
ory. Rigorous estimates for the growth of superficially convergent Green functions, 
and hence the p(k), can sometimes be obtained using constructive field theory, at 
least as bounds for the radius [26]. We emphasize that such results can be turned 
by our methods into similarly rigorous results for superficially divergent Green 
functions. 

A more complete discussion, dedicated to the renormalizable quantum field the- 
ories in four dimensions, will be given elsewhere. 

Appendix A. Other theories combinatorially with one primitive per 

LOOP ORDER 

For each of the following systems the solution series have coefficients satisfying 
the universal law. In the mixed 3 , </> 4 case there is one primitive per vertex per 
loop order. Unfortunately the full power of symmetry factors is not available in 
this simple combinatorial set-up, leading to the different variants. 
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Theory System Radius 





V 


— 1 + Aifi 




Smallest positive real root 






= xX 1 2 F|(l + Fi) 




of 3581577a; 4 - 4443984x 3 + 




x 2 






2332368t 2 — 5391 36r 4- 




2 




32768. Numerically 




F 2 


= 1 + F 2 X 2 




0.09061681898407704... 




Xi 


= i + x 2 xMf 1 + ' 6 - 


xX x F% 






1 1 


= xX x Fi{l + F 1 ) 




Root of a degree 10 
polynomial. Numerically 




^2 


= xF 2 + x 2 X x Fl 




u.izyoooyzzyouiy i ou . . . 






= 1 + F 2 X 2 






<f) A variant 




= 1 + x 2 X 1 3 F 2 i F 1 + xXiFi 
rest as in previous case 


Root of a degree 9 poly- 
nomial. Numerically 
0.13856076790723086 . . . 






=l+xX 1 X 2 F 3 z +xX'(F 3 i 








+ 2xF 2 X x 










+ X 2 X 2 F 3 (F a F b + 


F a + F b 








+ F a xXjFl + 


F b xX 2 F 2 








+ x 2 X\X 2 Fl - 


f xXfF* 








+ xX 2 F%) 










=1 + X 2 (F a F b + F a H 
+ F a .xX 1 2 F|+ 


-F b 








F b xX 2 F 2 + x 2 X 2 X 2 F£) 




mixed <fi 3 , cf> 4 




+ X*F 3 (F a F b + F a 
+ F a xX 2 Fi + 
+ x 2 X 2 X 2 Fl 


+ F b 

F b xX 2 F 2 


0.02145.. . 






+ xXfFl + xX 2 F 2 ) 





3 , 
+ ^X 2 F 2 

F a =x 2 X 2 F$ + xX 2 F 2 F a 
F b =x 2 XfFl + xX 2 F^F b 
X 3 =xX x F 2 + xF 3 + x 2 X 2 Fl 
F 3 =1 + X 3 F 3 
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